arXiv:l507.00910v2 [gr-qc] 26 Nov 2015 


On Stationary Axially Symmetric Solntions in Brans-Dicke Theory 


Pinar Kirezli^’^’* and Ozgiir Delice^’^ 

^Department of Physics, Faculty of Arts and Sciences, 

Namik Kemal University, Tekirdag 59030, Turkey 
^Department of Physics, Faculty of Arts and Sciences, Marmara University, Istanbul 34722, Turkey 

(Dated: November 30, 2015) 

Stationary, axially symmetric Brans-Dicke-Maxwell solutions are reexamined in the framework of 
the Brans-Dicke (BD) theory. We see that, employing a particular parametrization of the standard 
axially symmetric metric simplifies the procedure of obtaining the Ernst equations for axially sym¬ 
metric electrovacuum space-times for this theory. This analysis also permits us to construct a two 
parameter extension in both Jordan and Einstein frames of an old solution generating technique 
frequently used to construct axially symmetric solutions for BD theory from a seed solution of gen¬ 
eral relativity. As applications of this technique, several known and new solutions are constructed 
including a general axially symmetric BD-Maxwell solution of Plebanski-Demianski with vanishing 
cosmological constant, i.e. the Kinnersley solution and general magnetized Kerr-Newman-type so¬ 
lutions. Some physical properties and the circular motion of test particles for a particular subclass 
of Kinnersley solution, i.e., a Kerr-Newman-NUT-type solution for BD theory, are also investigated 
in some detail. 

PACS numbers: 04.50.Kd,04.20.Jb 


* pkirezli@nku.edu.tr 
1 ozgur.delice@marmara.edu.tr 



2 


I. INTRODUCTION 


Brans-Dicke (BD) scalar-tensor theory [1] is the most studied alternative theory of gravity generalizing general 
relativity (GR) in a consistent way by introducing a scalar field replacing Newton’s gravitational constant. Being one 
of the most straightforward extension of GR, this theory draws a lot of interest and helps to test various aspects of GR. 
The peculiar differences were noted in the exact solutions of this theory compared to the similar solutions of GR. For 
example, static spherically symmetric vacuum solutions of BD theory [2], unlike GR, do not describe asymptotically 
flat black holes [3] resulting from gravitational collapse, unless the scalar field becomes a constant. This is due to the 
fact that corresponding solutions in BD theory cannot meet three important criteria, namely, asymptotic flatness, 
regularity at the horizon and the weak energy condition, simultaneously. After the discovery of spherically and axially 
symmetric exact solutions, it was realized that there are ranges of parameters in which these solutions describe black 
holes [4-8] where the first two criteria are met but the third one is not. However, these ranges are at the nonphysical 
negative values of the BD parameter, w, where the kinetic term becomes negative. Hence, it is concluded that the 
black holes of BD theory obeying these three conditions must be the same as the black holes of GR, requiring a 
constant BD scalar. Although they are the same, their perturbations can behave differently [9] owing to the fact that 
they are the solutions of different theories. Thus, the differences between perturbations of GR and BD theories for 
the same black hole solution might be another tool [9, 10] to test BD theory against GR. 

Obtaining exact solutions of any theory is important for several reasons, such as, for comparisons with observational 
results or for obtaining the outcomes of the theory under consideration. One important class of these solutions with 
a great physical importance are the family of stationary, axially symmetric solutions, since the gravitational field of 
compact celestial objects such as stars, galaxies and black holes can be represented by such solutions. Investigating 
their properties may have some important astrophysical effects such as the existence and stability of orbits of planets 
and stars around such objects. In order to discuss these important properties of a gravitation theory, exact or 
approximate solutions must be obtained. For the BD theory we are considering, the field equations are more complex 
then those of GR, due to the presence of the extra scalar field. Even in GR, several solution generating techniques are 
developed due to the complexity of the field equations. Thus, in order to obtain exact solutions of BD theory, similar 
techniques were used. One of the methods frequently employed is making use of the Ernst equations [11, 12] derived 
from BD held equations [13-15]. Another method considers conformal transformation properties of scalar-tensor 
theories and using the equivalence of vacuum or Maxwell BD theory with Einstein-Maxwell scalar solutions [16, 17]. 

Besides these methods, following the earlier works on this subject [18-22], an identihcation technique is presented 
by Nayak and Tiwari [23] for axially symmetric vacuum solutions and generalized for Maxwell vacuum helds by 
Sing and Rai [24]. The main idea of this method is to employ a standard stationary, axially symmetric spacetime 
metric and the held variables sharing the symmetries of this metric and compare the held equations of both theories 
to hnd some transformations to reduce the held equations of BD theory to those of GR. Once the transformations 
linking general axially symmetric solutions of GR theory to BD theory are found, one can ehortlessly obtain new BD 
solutions for a given axially symmetric solution of GR. Thus, this method, which we call the Tiwari-Nayak-Singh-Rai 
(TNSR) method, is the most direct one. Using this method, Kerr-Newman [7] and black dihole type [25] solutions 
were constructed. On the other hand, one caveat for the present form of the TNSR method is that it contains 
only one extra parameter, the BD parameter, w, in addition to the parameters of the original GR solution. However, 
various vacuum or Maxwell vacuum solutions such as the spherically symmetric BD vacuum solutions [2], the Kerr-like 
solution presented in [19], and cylindrically symmetric static vacuum or Einstein-Maxwell solutions [26] all contain an 
extra arbitrary parameter in addition to w. Moreover, a recent work [27] presenting cylindrically symmetric Einstein- 
Rosen-type vacuum solutions also employed a similar technique involving two parameters. This fact implies that the 
TNSR method is not in its most general form. 

In this paper, our main aim is to investigate general stationary axially symmetric solutions of BD-Maxwell theory 
by Hrst obtaining a more general extension of TNSR method with one extra parameter that can be regarded as a 
measure of departure from GR. In order to do this, we first analyze the field equations for the general stationary axially 
symmetric spacetime given in the cylindrical coordinates in Sec. H. By employing a nonstandard reparametrization 
for the metric functions of a standard axially symmetric metric, we are able to obtain the Ernst equations from the 
field equations more easily. Analyzing this equation and remaining field equations together permit us to construct a 
two parameter extended version of the TNSR method in Sec. HI, in both Jordan and Einstein frames. Using this 
extended method we will first give some simple examples to show how this method works. Then we will obtain a 
BD version of the general axially symmetric Einstein-Maxwell type D solution, known as the Plebanski-Demianski 
solution in the case of the vanishing cosmological constant [28-30] in Sec. IV. We will also investigate some physical 
properties and the geodesics of a particular subclass of these solutions, namely, the BD version of a Kerr-Newman- 
NUT type solution, where the detailed calculations of the geodesic motion will be given in Appendix A. Sec. V is 
devoted to the discussion of ways of obtaining magnetized solutions for BD theory from a seed solution. Thanks to 
the technical convenience provided by this method, we have obtained the general magnetized BD Kerr-Newman-type 
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solution with a little effort using the corresponding GR solution and the extended TNSR method. As a last section 
in the paper, we will discuss the GR limit of these solutions in detail. This paper ends with a brief discussion. 


II. FIELD EQUATIONS AND SOME SOLUTION GENERATING TECHNIQUES 

A. Spacetime and field equations 

In this paper we consider Brans-Dicke scalar-tensor theory with a Maxwell field, described by the action in the 
Jordan frame as 


SjBD = J . ( 1 ) 

The field equations of this action can be expressed as 

+ ^ , ( 2 ) 

(2w -f 3) = T = 0. (3) 

Here is the energy-momentum tensor of the Maxwell field given by 

V = 2 , (4) 

and T = is its trace which is vanishing in four dimensions. 

A general stationary axially symmetric spacetime can be represented with a metric in cylindrical coordinates in the 
canonical form as 


ds^ = -e^^{dt + Adipf + + dz^) + e-^^W^dip"^, (5) 

where the metric functions U,K,A and W are the functions of p and z. This metric admits two Killing vectors dt 
and d,p. This metric can also be expressed in terms of other relevant coordinate systems such as Boyer-Lindquist 
coordinates by some suitable coordinate transformations. We also consider a Maxwell field sharing the symmetries of 
spacetime given by the following potential one-form: 

A = Ao{p, z) dt + Asip, z) dij). (6) 


B. Ernst equation for Brans-Dicke-Maxwell theory 


Now we discuss how to reduce BD field equations into the Ernst equation for axially symmetric spacetimes. In 
order to do this, following [31, 32], it is appropriate to use a different but equivalent metric Ansatz given by 

n n 

ds = —oe 2 (dt J- Ad(p( -t- oe ^ dp -t- —^ i,dp T dz ), (7) 

Va 

which simplifies the forthcoming equations considerably. The relation between the metric functions of (5) and (7) are 
given by 


Q 1 HI 

U = -1— In a, K = - \- ly -\— In a, W = a. 

4 2 4 4 


( 8 ) 


Note that although W and a are the same, we will keep both symbols and use whichever one is relevant, in order to 
remind ourselves which metric Ansatze we are using. 
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In the presence of a Maxwell field of the form given in ( 6 ), BD-Maxwell field equations can be written in a compact 
form as 


= 0 , 


iv. + 2e^/'^[A - V^oVAg] = 0, 

iv. {a(t)V^ +AV. (a<j>e^VA^ +a<j>e^{VA)^ 

+ 2e9 A'^ (vAoY - (vAsY - 2e-9 (vAq)" = 0, 
(vm) - ^acj) ( 

^^(2a; + 3)V (aV^ = 0, 


Vfl) V - 2a(/)VV- 


(V 0)2 „ 

' — uja — 7 — = 0, 


V. 

V. 


g-n/2y^o + (vAs - AVAoj 

gn /2 


= 0 , 


= 0 . 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


The above equations (9)-(12) are obtained from BD field equations (2), Eq. (13) is the scalar field equation (3) and 
(14) and (15) are the nontrivial Maxwell equations which do not involve the scalar field (f). Here the differential 
operators are the operators in the fiat Minkowski spacetime, i.e., V = dpC^ + such that = 9^+9^+ 9^. 

Note that in Eqs. (9)-(ll) the functions a and (f) are present only as products acj). Following [11, 12, 32] we can put 
Eqs. (10),(11),(14) and (15) into two complex equations given by 


(Re £ + |$n-V.(dV£:) = (Ve + 2$*V$).V£, 
a 

(Re e + |$P)-V.(dV$) = (Ve + 2$*V$).V$, 
a 

where the new function a is given by 


a = (jja = 4>W, 


(16a) 

(16b) 


(17) 


satisfying Eq. (9), i. e., = 0. The complex potential $ and the complex function £ are defined as 


e = / - + ih, 


(18) 


$ = Ho + fi3. (19) 

Here the scalar function / is given by 

/ = ( 20 ) 

and new vector potentials A 3 and h are defined by 

epXVA 3 = e^/^{VA 3 -AVAo), ( 21 a) 

4 X V/i = e^dVM - 24 X Im($*V$), (21b) 

where belongs to the set of basis vectors of three-dimensional flat space {e?, e^, e^}. 

Equations (16a) and (16b) are actually in the same form as the corresponding equations in GR [33]. The real and 
imaginary parts of the equation (16a) are equivalent to the equations ( 10 ) and ( 11 ) whereas the real and imaginary 
parts of Eq. (16b) are equivalent to Eqs. (14) and (15). 

The field equation (9), namely, V^(q;(/)) = V^d = 0, permits us to employ the canonical value for d, such that we 
can set a = p. This brings Eqs (16a) and (16b) into 

{s + s* + |$| 2 ) v 2 £ = 2 (v£ -b 2 $*V$) .V£, 

(£ -b £* -b | 4 >P) = 2 (v£ -b 2 $*V$) .V$, 


( 22 a) 

( 22 b) 
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where now is the Laplacian operator in the three-dimensional flat cylindrical coordinates. These equations are 
exactly in the same form with the Ernst equation [11, 12] for GR in the presence of a Maxwell field. Thus, the 
choice of the metric form (7) simplifies the calculations in obtaining Ernst equations for BD theory considerably. 
Actually, the above Ernst equations for BD theory were already derived before [13-15] in the Einstein frame, and 
some exact solutions were obtained by integrating these equations. Note that the above Ernst equations of BD theory 
in the Jordan frame can also be put in the the other useful forms of the Ernst equations known in the literature with 
exactly the same transformations. However, we will not pursue this path further in this paper. Note also that, for 
the following discussion, we will retain the a term in order to cover some exact solutions where a ^ p. 

For given Ernst potentials solving these two equations, the remaining field equations (9),(12) and (13) must be 
integrated to obtain the full solution. The fact that the form of the BD Ernst equations given above are exactly the 
same equations with the Ernst equations of GR enables us to make an important observation that any solutions of the 
Ernst equations of GR are also the solutions of BD theory formulated in the manner presented above. This requires 
the identifications 


a = a(j) = aE, = ^e, A = Ae, A = Ae, (23) 

where subscript E represents any solution of the Ernst equation in GR. The above relations are due to the unconven¬ 
tional parametrization ( 8 ) of the stationary axially symmetric metric (7). Note that there is an arbitrariness in the 
metric function a and the scalar field 4> for the choice a = acj), since any choice of a and </> satisfying Eq. (9) is not 
guaranteed to solve the equation (13). Hence, among all the possible choices of a and (f) satisfying (9), only the ones 
satisfying Eq. (13) are the solutions of the BD field equations. Thus, although it is a powerful method, obtaining 
solutions from Ernst equations still requires some effort to solve the field equation (13) and integrate Eq. (12) to find 
v. Hereafter, we will employ a more direct method, yielding the solutions readily from known GR solutions, up to 
possible coordinate transformation. Actually, as we have discussed in the Introduction, there is such a method in the 
literature, given in [23, 24], which successfully generated many BD(-Maxwell) solutions from the known solutions of 
GR theory. In the next subsection, a more general form of this method, containing one extra parameter apart from 
the BD parameter w, will be presented, with the help of the Ernst equations obtained in the above analysis. 


III. EXTENDED TIWARI-NAYAK-RAI-SINGH METHOD 

It is usually difficult to directly solve the field equations of Einstein’s or Brans-Dicke theory of gravity to obtain 
exact solutions. Due to this difficulty, some solution generating techniques are developed. Here we extend one of them, 
which enables one to generate an axially symmetric stationary or static solution of Brans-Dicke(-Maxwell) theory from 
a known solution of Einstein(-Maxwell) gravity with the same symmetry properties. The essence of the method is 
the following. One starts with general BD(-Maxwell) field equations of a metric of the form (5), where the scalar field 
(j) and the Maxwell field share the symmetries of the metric, with metric and field variables K, U, A, W, A, (j). Then 
by investigating the transformations on these metric functions such that the field equations exactly reduce to field 
equations of Einstein(-Maxwell) theory for the metric and field variables of GR with Ke,Ue, Ae,We, Ae- 

An immediate observation of the field equations (9)-(I5), as we have discussed above, is that in Eqs. (10) and (II) 
or equivalently in (I6a) and (I6b) the metric function W{= a) and the scalar field (j) enters the expressions only as 
products a = acj), owing to the form of the metric Ansatz (7). Since, in the above choices we have considered, the form 
of the Ernst equations are similar for GR and BD theories, any solutions of this equation for GR can be considered as 
solutions of this equation in BD theory provided that the product a(j) is equal to the corresponding metric function 
for Einstein-Maxwell theory ue, together with the remaining metric functions and Maxwell fields except v. Note that 
for this choice at/) = ue, it is easy to show that the requirement that a and (j) solve the remaining field equations (9) 
and (13) is met by the Ansatze a = </> = with k being a constant parameter, measuring the departure from 

GR. Actually, the choice is well motivated, since this Ansatz is a general solution for the canonical value ue = p with 
the choice </> = /(p), a = pjf{p). For these choices Eq. (13) becomes 

y) = 0 ’ ( 24 ) 

whose general solution is 

f{p)=P^\ (25) 

and we have chosen the integration constant cq = 1 —fc similar to some previous solutions [26, 27] exist in the literature. 
We will not consider the possibility of more general choices of a and (j) in this paper. 

By a careful analysis we have obtained the following theorem. 



6 


Theorem 1. From any solution of the Einstein or the Einstein-Maxwell theory with the metric of the form (7) labelled 
by a subscript E and a possible Maxwell field for the Einstein-Maxwell case, i. e., Ae = Aq{p, z)dt + A^{p, z)d(f>, the 
corresponding BD or BD-Maxwell solution can be obtained by the following transformations: 


ex — ttg. 

(f = a^E^ , Ft = Qe, 

(26a) 

A = Ae, 

A = Ae, 

(26b) 


2c^ — (2c(j -|- 3)/p 


V = ve + 

4 

(26c) 


To prove this theorem we only need to prove the relation (26c), which can be easily obtained by applying the 
transformations (26a) and (26b) to BD-Maxwell field equations above and solving for ly. Thus, for this form (7) of 
the axially symmetric stationary spacetime, the transformation equations relating GR solutions to BD ones turn out 
to be quite simple. The advantage of this form of the metric is that the Ernst equations of BD and GR theories 
are obtained straightforwardly in these coordinates and the related developments for the Ernst equation are easily 
applicable to BD theory in this parametrization. However, the form of the metric (7) is rather unfamiliar. Hence, 
in order to compare our results with the existing solutions in the literature and also to apply this method to obtain 
new solutions, we need to express the above theorem in terms of the metric (5). By considering the relations (8) the 
above theorem becomes as follows. 


Theorem 2. From any solution of the Einstein or the Einstein-Maxwell theory with the metric of the form (5) labeled 
by a subscript E and a possible Maxwell field for the Einstein-Maxwell case, i. e., Ae = Aq{p, z)dt A'i{p, z)d(f>, the 
corresponding BD or BD-Maxwell solution can be obtained by the following transformations: 


= A = Ae, A = Ae, 


U = Ue 



K = Ke + 


2cc — 1 — (2cu -f 3)/c 
4 


In <j). 


(27) 


The only possible technical difficulty in using the above method is to put the GR solutions into the particular form 
(5), since the stationary axisymmetric solutions can be expressed in various coordinate systems. For convenience, let 
us present the corresponding line element in the Jordan frame explicitly as 

ds'^jBD = e^^’^idt + AEdiff + SKE-UE)^{i-k){u+3./2)i^^^2 ^ ^ (28) 


Now we can compare this result with the results present in the literature. Note that the above theorem is known in 
the literature in a more restricted form, in which in the works [23, 24] the case K = Ke is considered, which relates 
the parameter k with the BD parameter, oj as k = (2a; — l)/(2a; -I- 3). The restricted form with K = Ke is used 
to construct Kerr-Newman -type solutions [7] in BD theory from GR and also Bonnor-type black dihole solution 
is discussed in [25]. Note that a similar construction was used to generate Einstein-Rosen-type gravitational wave 
solutions in BD theory [27] for the special case W = p. 


A. Extended TNSR method in the Einstein frame 


It is well known that a conformal transformation of the form 


— 4 ^ 



\a(j) 


brings Jordan frame BD action into the Einstein frame as 


Seed = J d^x^g (^R - . 


(29) 


(30) 


Note that we use the unit system where SttG = 1. One important property of the above transformation is that the 
Maxwell Lagrangian is invariant under the conformal transformations in four dimensions. Hence, unlike in the case 
of other matter fields, the scalar field is not coupled with the Maxwell field and the Einstein frame of BD theory 
for the Maxwell field is equivalent to Einstein-Maxwell-scalar solutions. Now, let us investigate the results of the 
previous section in the Einstein frame. Recall that for a given axially symmetric solution of Einstein-Maxwell theory 
has a corresponding BD-Maxwell solution given in Theorem (2). After applying the conformal transformation (29) to 
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the metric (28) and the scalar field, the corresponding solution of BD-Maxwell field equations in the Einstein frame 
becomes 


+ Asdipf + ^- 2 Ue ^2 ^^2 

, (l-fc)(2w + 3), 

^ = 2 -- '-\nWE- 


Clearly, we have the following theorem. 


(31) 

(32) 


Theorem 3. From any solution of the Einstein or the Einstein-Maxwell theory with the metric of the form (5) labeled 
by a subscript E, and a possible Maxwell field for the Einstein-Maxwell case, i. e., Ae = Aq(p, z)dt + A^{p, z)d(j), a 
corresponding BD or BD-Maxwell solution in the Einstein frame (Einstein-Maxwell-scalar solution) can be obtained 
by the following transformations: 

W = We, + x^We, A = Ae, A = Ae, (33) 

U = Ue, K = Ke + 

Note that a solution generating technique was presented for Einstein-Maxwell-scalar (EMS) theory in [22] where 
in that work it was shown that, for a given Einstein-Maxwell solution, the corresponding EMS solution can be found 
by solving a set of differential equations involving only K and if in our notation. 


B. Examples in fordan Frame 

1. Flat spacetime 

Let us apply the above procedure to flat Minkowski spacetime, to demonstrate the method and clarify some issues 
in the previous restricted form of this technique. If we consider Minkowski spacetime in Cartesian coordinates, we see 
that since We = 1, the resulting solution is also the same flatspacetime in Minkowski coordinates, since </> = 
becomes constant. In order to obtain a nontrivial scalar field, we have to use coordinate systems where We is not a 


constant, such as the flat spacetime given in spherical coordinates: 

ds^ = -df + dr^ + r^{de^ + sin2 OdgA). (34) 

Note that here r denotes spherical radial coordinate. In order to bring this metric into the form (5), we can consider 
the transformation r = e^, which brings the metric into the form as 

ds"^ = —df^ -\- e^^{dR^ -\- dO^ -\- sin^ Odip^). (35) 

Here, by comparing this with (5), we see that it resembles the form with the identification R = p, 9 = z, and 

WE = e^sm9, Ue = 0, Ke = R, Ae = Ae = 0. (36) 

Now applying the procedure in Theorem 2, we obtain 

ds2 = {rsiTi9)^-^[- df + (rsin6»)(i-'=)'(‘^+3/2)(dr2 + {r sin 9)^de)^], (37) 

<) = (rsin9)^-^ (38) 


Clearly, the application of this method breaks the spherical symmetry of the BD empty spacetime, as discussed in 
[25]. Namely, even if one starts with a spherically symmetric static solution as a seed, the resulting metric will be 
also static but axially symmetric. Thus, this method is not appropriate for obtaining the spherically symmetric BD 
solutions and one may use the methods given, for example, in [17] to obtain spherically symmetric solutions. The 
reason is that the g^^ and g^p^p components are changed after the application of this method. 

One can also consider the cylindrical flat spacetime given by 

ds^ = —dt^ -\- dp^ -\- dz^ -\- p^dip^, (39) 

as the seed solution with Ke = Ue = 0, = p. After applying the algorithm, the resulting flat space solution of 

BD theory becomes 

ds^ = [-de + p(l-0=(-+3/2) (^^2 ^ ^^2) ^ ^2^^2j ^ (40) 

Note that one can also obtain this solution from (37) by usual transformation relations between spherical and cylin¬ 
drical coordinates. Clearly, unlike spherical symmetry, the cylindrical symmetry is preserved by this method. 



2. Bonnor type dipole solution 


As we have discussed, this method cannot be used to construct spherically symmetric solutions even if one starts with 
a spherical seed solution. The resulting solution becomes axially symmetric. Thus, we cannot obtain a Schwarzschild- 
type solution using this algorithm. However, for static axially symmetric solutions we can use this method to obtain 
new solutions with the same symmetry. Using this fact, we now generalize a previous BD solution constructed using 
the restricted form of the TNSR algorithm to its more general form. The solution we choose as the seed solution 
is a static axially symmetric solution of Bonnor [34] representing a magnetic dipole in GR. The corresponding BD 
solution of the magnetic dipole presented by Bonnor as a solution of Einstein-Maxwell theory was presented in [25] 
using the TNSR algorithm. Since the original method is restricted, we now present the more general version of this 
solution using the extended version we have discussed. The GR solution reads 


ds^ = ^1 

Aq = 0 , 


2Mr\ 


— dv + 




[A + (M2 + a2) sin^ 9 " 


dr"^ 


A sin^ 6 


(1-^) 


:d(p^, 


A,= 


2aMr sin^ 6 
A + a2 sin^ 9 ’ 


A = - 2Mr + 


Yi = cos^ 9. 


(41) 

(42) 


In order to apply the algorithm, we need to bring the metric into the form of (5). To do this we need to apply a 
transformation r = + M + (M^ + a^)e“^, bringing the {r,9) sector of the metric into the form dr'^/A + d9^ = 

dR^ + d9'^. Hence by comparing the obtained metric with (5), we can read the metric functions of Einstein-Maxwell 
theory as follows: 


g2C/^ 





Ae = 0 , 


1U| = A sin^ 9 


[a -|- (M2 -I- a2) sin^ 9 ^ 



4 


(43) 


By applying the extended TNSR algorithm and transformation back to the Boyer-Lindquist coordinates, the BD 
version of the Bonnor’s solution becomes 


ds^ 

4> 


A sind 


fe-i 


i-fc 


1 - 


2Mr\‘ 


— dt^ -\- 


A sin 


in6») 


{l-kY{uj+Z/2) 


2 


[A 4- (M2 -I- a2) sin^ 9 


Asind) , A = A^dLp. 



A sin^ 9 






, (44) 


(45) 


Note that this solution reduces to one given in [25] for the restricted case given by fc = {2uj — l)/{2uj -\- 3). 


IV. A GENERAL AXIALLY SYMMETRIC SOLUTION: BRANS-DICKE-KINNERSLEY SOLUTION 

A. Brans-Dicke-Kinnersley solution 


The general Einstein-Maxwell type D solution, known in the literature as the Plebanski-Demianski solution [28], 
first presented by Debever [29] is given by the metric [33] 


ds^ = (1 - pq) 


-2 


X [dt + q^da)"^ — Y{dt— p^daY 

p2 _|_ q2 




Y 


X = {—g^ + 7 — A/6) -I- 2lp — ep^ + 2mp^ — {e^ + y + A/&) 
Y = (e^ -b 7 - A/ 6 ) - 2mq + eq^ - 2lq^ + - 7 + A/ 6 ) q^. 


(46) 

(47) 

(48) 


Here m is mass, I is the Newman-Unti-Tamburino (NUT) parameter, 7 and e are related to the angular momentum 
per unit mass a, and the acceleration b and e and g are the electrical and magnetic charges. Now we want to find 
the BD version of this solution. Thus, we need to cast the metric to the form given in (5). To do this, since the 
method do not work in the presence of a cosmological constant, we set A = 0 [30]. We also factorize the dt, da sector 
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of the metric and apply the following coordinate transformations dp = \/Xd6, dq = y/Ydr, a = (p, then the metric 
components of GR for this solution in the form (5) becomes 


2 KE = \n 
2 C/b = In 


Y - X 


(1 -M)4_ 

Y-X 


Ae = 
We = 


q^X +p^Y 
X-Y ’ 
y/XY 
{l-pqf 


(49) 


[{l-pqY + g2)J 

By applying the method given in Theorem 2, and transforming back to the p, q coordinates, the BD version of the 
Plebanski-Demianski solution with a vanishing cosmological constant is found to be of the form 

I k—1 


ds^ = 


yfXY 




(1 -m) 


X {dt + q^da) —Y{dt — p^da) ‘ 
p^ + q^ 


fdp^ dq^\ 

y/XY 

(l-fc)2(aj+|) 'I 

1 rh — 

y/TY 

(^ + —) 

_{^-pqf _ 


_{^-pqf _ 


l-k 


(50) 

(51) 


As to the best of our knowledge, this solution for BD theory is new. Note that, setting the acceleration parameter 
to zero and employing Boyer-Lindquist-type coordinates, one can recover a BD Kerr-Newman-NUT-type solution. 
Note alsothat, due to Hawking’s theorem [3], these solutions do not describe asymptotically flat black holes originated 
from the gravitational collapse of stars, unless the scalar field is a constant. However, these solutions may represent 
the exterior fields of compact bodies whose interior is characterized by the parameters of this solution. Hence, it is 
still important to investigate certain properties of this solution. 


B. Kerr-Newman-Taub-NUT solution in BD theory 


Here we obtain the BD version of general Kerr-Newman-Taub-NUT (KNTN) solution using the extended method. 
Note that the solution can also be obtained from the Plebanski-Damianski-type solution discussed above by a relevant 
limiting procedure. In GR, this solution has the metric and Maxwell field as follows: 


ds^ = — 


A — sin^ d 


{dt + Ad(j))^ 


A sin^ 6 


(r^ -I- -I- — abY 

A — sin^ 9 


d(j) 


dr^ 

~A 


+ d9^ 


Ae = - bd<j)), 


where the metric functions are given by 

A _ 2 r> I 2 2 I 2 

A = r — 2m r + a — n + q , 


Y = r'^ + {na cos 0)^ 


(52) 

(53) 

(54) 


A = 


a{r'^ + a'^ + ri^) sin^ 9 — bA 


By transforming the radial coordinate as r = + m+ (m^ — a?' -\- 


A — a? sin^ ( 
— . 


h = a sin"^ 6 — 2n cos 0. 


■ n* — (?^)/4 e 

?2 I 


-R 


and using the transformation given 


in [35], the {r,9), sector of the metric becomes dr^ / A -f = dR^ + d9'^ such that it resembles the form of the metric 
(5) with the identification R = r, 9 = z. Then, one can apply the above theorem and find the corresponding solutions 
in BD theory easily in these coordinates. Transforming back to the original coordinates, the KNTN solution in BD 
theory can be expressed in Boyer-Lindquist-type coordinates as follows: 


ds^ = (-x/Asind) 


k-l 


A — sin^ 


{dt + A d(j)y 


A sin^ 9 


(r^ -I- -I- — abY ^ 

-d4>^ 


A — a? sim 


Asind) + 


(55) 


A = Ae. (56) 

This solution is also presented in [8] by using a sigma-model analysis with having a different constant parameter a 
where the relations of this parameter and the parameter k are given by a = (1 — k){2oj + 3)/4. The Ricci scalar of 
this solution is given by 


R= {k — Ifuj 


2 AC0S^9 + {M-rYAxYd f fe-3-i(fc-l)[l-2c;+fc(2c^+3)] 


A sin 


ind)' 


(57) 
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which is helpful to better understand spacetime singularity structure and it vanishes for the limit fc —>■ 1. Clearly, 
there is a ringlike singularity at the location described by S = 0. The last term in the parentheses vanishes if its 
exponent becomes zero. Note that the parameter range of this spacetime is investigated [8] for the existence of a black 
hole type solution with a regular horizon by considering some regularity and existence conditions and concluded that 
for w < —3/2, a± = (1 ± a/ 1 + 4|2a; + 3|)/4 spacetime may represent a black hole with a regular horizon. In our 
notation the latter condition is = 1 — (1 + a/ 1 + 4|2w + 3|)/(2a; + 3). Note that the symmetry axis is non-null 
for a_ but null for a_|_. In these black hole-type solutions, there is a caveat: the scalar field do not satisfy the weak 
energy condition. Hence, these black holes do not satisfy Hawking’s criteria [3]. Moreover, this black hole has a 
vanishing surface gravity and they are ergocold since they have zero Hawking temperature [7, 8]. 


C. Circular geodesics of the BD-KNTN Spacetime 

Most of the previous studies of axially symmetric spacetimes in BD theory are focused on the obtaining the solutions 
and analyzing the solutions for whether these solutions represent black holes. Another property of these solutions 
with possible astrophysical importance is the motion of the test particles. There are few works on this issue [36]. In 
this paper we have also investigated the circular geodesics of BD-KNTN solution. However, since the calculations 
and equations are rather complicated and lengthy, we present them in the Appendix. Here below we summarize our 
general resuls. 

We have first calculated the general geodesics equations in the case of equatorial motion. Then we have considered 
null and timelike geodesics. For photons, we have reduced general radial geodesics equations to integrals depending 
on the radial coordinate in the special case L = aE where E is the energy and L is the angular momentum of the test 
particles or photons. For general case, we have been able to obtain the equation of the photon sphere and realized that 
it is the same as the corresponding GR case. For timelike particles, we have obtained an effective potential equation 
for radial motion. We have analyzed this equation for the existence of the inner most stable circular orbit (ISCO). 
After long and complicated calculations, we are able to obtain the equation determining ISCO. This equation has 
been analyzed for the special cases such as the extremal case and the BD-Kerr case, as well. These equations agree 
with the corresponding limiting cases in GR. 


V. MAGNETIZED SOLUTIONS IN BRANS-DICKE THEORY 


Obtaining a magnetized solution from a given seed solution using Ehlers-Harrison-type transformations has a long 
history [37]. Using these transformations one can generate a magnetized solution from a given Einstein-Maxwell 
solution. For example for a given axisymmetric solution with property gi^p = 0, the following transformations [38] 
bring a given solution of EMS solutions (i.e., the BD-Maxwell vacuum solutions in the Einstein frame) into the 
magnetized one: 


= A2 




4' =_fU 

^ KB 


9^^ ='^9cpcp,a'= a, 


1 -I- 


B- ^ = 


1 -I- -j^BAp 


+ p^B 9p(p- 


One can also transform the obtained solutions to the Jordan frame by using (29) to obtain magnetized BD solutions 
as well. This method is used for BD theory [25] in connection with the previous version of the TNSR method [23, 24] 
to obtain a magnetized black dihole-type solution. However, if a magnetized solution is already presented in the 
Einstein-Maxwell theory, we can use the extended TNSR method to obtain the corresponding magnetized solution in 
BD theory for Jordan or Einstein frames considering Theorem 2 or 3. 

Using the extended TNSR method discussed above, we can now present some new magnetized BD solutions such 
as the magnetized Kerr-Newman-type solution or the magnetized Bonnor dihole solution. 


A. Magnetized Kerr-Newman solution 


The magnetized Kerr-Newman solution is first discussed in [37, 39] but since the Harrison transformations for 
stationary metrics are more complicated and cumbersome then the static ones, the full exact solution is first presented 
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recently by [40]. In their notation, the solution reads 

' dr'^ 


ds = H 


-fdt^ + 


de^ 


Ae = ^o + $3 {dip - Qdt) 


where 


E sin^ 9 

HR^ 


{dp — Vtdt)^ , 


(58) 

(59) 


+ a? cos^ 9 , A = (r^ + a^) - 2mr + 

f=^, J: = {r^+a^f-a^Asin^9, (60) 

2-i 


and the (very lengthy) expressions of the other metric and field variables H,n{= ui in their paper), $o,4>3 can be 
found in [40]. 

To use the extended TNSR method, we need to bring the above metric into (5), and to do this we should regroup 
the t, p sector of the metric and consider the transformation dr/\fK = di?, then we can read off the metric functions 
of the GR solution as 


^2(7e _ 


= Hf- 


sin^ 9 

’ 




Ae = 


flE sin 9 


- En2 sin^d’ 


We = ■\/Asin0. (61) 


Then applying the extended TNSR method, the corresponding magnetized Kerr-Newman type solution for BD theory 
in the Jordan frame can be obtained as 


ds^ 




fdf + R^ 





^^{dp-ndtf^- ( 62 ) 




(63) 


Note that the Einstein frame solution can be obtained by considering only the terms in the curly bracket and replacing 
<j) with Ip as in (31) and (32). The physical properties of this solution will be treated elsewhere. 


B. Magnetized Bonnor dihole solution 


Now as a final application of the method discussed in this paper, we present a two parameter extended version of 
the BD solution [25] for Bonnor-type dihole solution [34] embedded in a Melvin magnetic universe presented in [41]. 
Starting from Einstein-Maxwell magnetized solution and applying the Theorem 3, we obtain the following solution: 


ds"^ 


A 

A 


A sin 9 


fc-i 




( \/A sin 0 j / j 2 

-dt^ + -^^ ^ + d9^ 

^ A 


[A -I- (M^ -|- a^) sin^ 0] 

V L 

_ \ 1 — k 

Asin^j , 

2Mra + 3/2 sin^ 6 


A^dp = - 


AE 


sin 9 dp, 


A + sin^ 9 + 2BMras\T? 9 + B"^/dsiv? 9 (r^ — -f Ao? sin^ 9 


A sin^ 9 
A 2 


dp^ 


(64) 

(65) 

( 66 ) 

(67) 


VI. GENERAL RELATIVISTIC LIMIT 

The GR limit of the solutions of the BD theory is not as trivial as it sounds. It was often claimed that as a; i—>■ oo 
BD solutions reduce to the corresponding GR solutions. However, several counterexamples of this claim are present 
in the literature. Eor example, for the special case of vanishing trace of the energy-momentum tensor, T, several 
papers [26, 42-46] are showed that this is not true in general. This fact is related in [44] with the conformal invariance 
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of BD theory for vanishing T. This may be also due to fact that the scalar field equation becomes sourceless and 
the integration constant of this scalar equation becomes independent of w since the remaining field equations cannot 
fix this constant in terms of w. A perfect example is a recent paper [46], in which BD-Maxwell solutions for higher 
dimensional static cylindrically symmetric spacetime were obtained. In these solutions, when the spacetime dimensions 
is greater than 4, the trace T does not vanish and this yields an equation relating BD parameter w and the extra 
integration constant of the solution and solving this constant for w, one recovers the corresponding GR solutions [47] 
in the w i—oo limit. However, in four dimensions, the field equations do not yield such an equation and for this case 
[26] the GR limit is achieved not by w i—>■ oo but a specific value of this integration constant where the BD scalar 
becomes a constant. Even for a nonvanishing T, some counterexamples were found [48] where the scalar field does 
not go to a constant sufficiently fast and therefore the GR limit is not recovered as w i—>■ oo for these solutions. 

For the cases where this limit do not arise as a result of field equations, some physical arguments can be put forward 
to fix this extra constant in terms of w. However, these arguments are not in general the result of encompassing physical 
requirements but the results of some special considerations. For example in [45] where some spherically symmetric 
BD solutions were discussed, the extra constant is fixed by demanding the matching of the solution to an internal 
source in the weak field approximation. Although the GR limit is obtained as w i—oo for this special case, one cannot 
ensure that this procedure also works for the full theory, especially in the presence of a source having a strong field. 
Thus, in order to investigate the GR limit of BD solutions, not only mathematical results but also additional physical 
arguments might be required, as in [45]. 

In order to avoid all these complexities in this paper, when we have discussed the GR limit, we have only considered 
the case where the scalar field becomes a constant. As it is clear from Eq. (26a), as fc —>■ 1, the stationary axially 
symmetric BD (-Maxwell) vacuum solutions reduce to GR for any finite uj irrespective of whether it is large or small, 
since </> becomes a constant. Note that, similar to the case of generic finite values of w where the BD theory is equivalent 
to GR with a minimally coupled scalar field, modulo potential singularities in transformation between frames, when 
UJ !->■ oo, it is known that the theory reduces to GR plus a minimally coupled (residual) scalar field [49] even if 
lim(^Hj.oo 4> = 4>o = constant. Hence, we need some consistency requirements on the k, since this requires lim(^Hj.oo k = 1. 

For example, the regularity of the spacetime can be ensured by choosing k such that lim^H_i.oo /c ~ 1 — with 

/3 an arbitrary constant. In the limit uj i— oo, for this value of fc, the metric function K does not approach to its 
corresponding GR value but becomes K = Ke + f InlRE. Hence, as w i—>■ oo, BD solutions discussed in this paper 
do not approach their corresponding GR solutions. Note also that although for k ^ 1 the GR limit is obtained 
since the BD scalar becomes a constant, the additional scalar degrees of freedom is always present at the level of the 
fluctuations, even for the constant background solutions. 


VII. CONCLUSIONS 

In this paper, we have investigated stationary axially symmetric vacuum solutions of the Brans-Dicke-Maxwell 
theory in both Jordan and Einstein frames. We first show that, employing a particular form of the standard metric 
simplifies the procedure of obtaining the Ernst equation from the BD field equations. By investigating the Ernst 
equation and also the remaining field equations, we are able to obtain a two parameter extension of a particular 
solution generating technique for BD theory. In order to show how this method works, we have constructed several 
known solutions and also some new solutions for BD theory, such as a Plebansky-Demiansky-type solution with a 
vanishing cosmological constant or magnetized Kerr-Newman-type solutions. We have also discussed the GR limit of 
these solutions in some detail. 

As a further physical application, we have also investigated the geodesics equations for Kerr-Newman-NUT type 
solution of BD theory in Appendix A, where we have focused on the circular null and timelike geodesics of this 
solution and successfully obtained the ISGO for this spacetime. A question arises whether this known solution 
generating technique generalized in this paper to two parameters can be applicable to other theories such as f(R)-type 
theories or EMS-type theories. Especially, EMS-type theories are quite interesting since a regular black hole type 
solution with scalar hair [50-52] exists and interest has recently been focused on the possibility of the existence of 
stationary-type hairy black hole solutions [15, 53-57]. 
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Appendix A: CIRCULAR GEODESICS OF THE BD-KNTN SPACETIME 

In order to determine the geodesic equations of BD-KNTN solution given in Sec. IV B, we first define the Lagrangian 
of the system as 


C 


1 / -2 • . .2 

2 \ 9ttt + gtipW + + gr 


(Al) 


where here the overdot denotes the differentiation with respect to the proper time t for timelike case and an affine 
parameter for the null case. The symmetries of this Lagrangian enable us to determine the following first integrals of 
the motion 


gtti + g^tV = -E, gt^i + gw^ = L, 


(A2) 


where the constants E and L are related with energy and angular momentum of the particle. Furthermore, using the 
metric itself, we have 

e = gui"^ + ‘^■gt^iip + g^^(p‘^ + grr'r'^ + geeO'^, (A3) 

where e = 1,0,—I for spacelike, lightlike and timelike geodesics, respectively. From (A2) we obtain 


t = 


ip = 


IT2 


e^^{L + EB) 

■ 


(A4a) 

(A4b) 


To determine the geodesics in the equatorial plane, we set 0 = | and 9 = 0, and also use (52), then the above 
equations become 


if = 

i = 


(A — a^)(L — aE) + aE{r'^ + n^) 

(r2 n2)A('=+i)/2 ’ 

[a(A — a^)(T — aE) + (r^ -|- n^){E{r‘^ + to?) — a{L — 2Ea))] . 


(A5a) 

(A5b) 


(j,2 _j_ ^2'j/\(k+l)/2 

Using these expressions in (A3), we obtain the radial equation determining the geodesics of BD-KNTN spacetime as 
Ae-2*^ + EB)^ + E^W^ -k elU^e^t^] 


= 


IU2 


(a^ — A)(L — aE)"^ + {r"^ + ri^) [U2(r2 -|- n^) -I- 2aE{aE — L) + eA(^+^)/2j 

(■^2 _|_ „2)2/^i=^[2c;-l-fc(2i^+3)] 


(A6) 


These last three equations are the general geodesics equations for equatorial motion. Now, let us consider null and 
timelike motion in some detail. These equations reduce to corresponding GR ones [58, 59] for fc i—>■ 1. 


1. Light-like geodesics (e = 0) 

Null geodesics are obtained by setting the parameter e = 0, then the radial equation becomes 

Ae-2^ [-e*^{L + EB)^ + E^W^ + e 1^2 


f2 = 


IU2 


(a^ — A)(L — aE)'^ -I- {r'^ + n^) [U2(r2 -|- ri^) -|- 2aE{aE — L)] 

~ (-^,2 _(_ ^2)2^i^[2(.,;-l-fe(2i^-|-3)] 

In order to distinguish the geodesics with different impact parameter D = L/E, we consider the following cases: 
(a) The special case L = aE: 


(A7) 
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For this case the impact parameter becomes D = a, which plays an important role to study the radial geodesics. 
Equations (A5a),(A5b and (A6) become 


aE 

^ ~ A('=+i)/2 

(a^ + r^ + v?)E 
* ~ A('=+i)/2 

(A8a) 

(A8b) 


(A8c) 


Here the positive sign represents outgoing photons whereas the negative sign represents ingoing photons. Using 
these equations we find that 


±ip = aj dr, (A9a) 

±t = J{a^+r^+ n2)AM('=-i)^(2<.+3)-8] 

± Ur = y A ^[ 2a ;- l - fc ( 2a ;+ 3 )]^^^ (^ g ^) 

These equations are more complicated than similar equations for GR case [58, 59]. For example, unlike GR 
case, the radial coordinate is not changing uniformly in BD theory. 

(b) The general case D = L/E: 

Using Eq. (A7), we obtain the equations corresponding the radius Tc of the circular photon orbit aX E = Ec 
and L = Lc as 


1 


, (a2-A)(aF;,-Le)2 


{ri+n^)E^ 


- 2aEc{Lc - aEc) 


{rl + n^) 

For the regions where A ^ 0, the derivative of the above equation becomes 


= 0 . 


(AlO) 


X L, [(re - m){rl + n") + re(a" - A)1 1 = 0. (All) 


Substituting Dc = ^ into the above equations we obtain 


, n 2^ (a^ — ^){a — Dc)^ „ ^ 

(r^ + n ) H - 2a{Dc - a) = 0, 

^ + n^) + - A)] = 0. 


(A12) 

(A13) 


Note that these equations are independent of the BD parameters k and w, meaning that the BD scalar field 
does not affect the geodesics of the photons. The impact parameter can be obtained from Eq. (A13) as 


Dc = a ± 


rdr"^ + 


{tc — + V?) + rc{a? — A) 


(A14) 


Inserting this equation into (A12) yields the equation of photon sphere [60, 61] for BD-KNTN spacetime, which turn 
out to be the same with the corresponding GR solution [58, 59]. 
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2. Timelike geodesics (e = — 1) 

For timelike geodesics, Eqs. (A5a) and (A5b) remain unchanged, but Eq. (A6) becomes 
(r^ + n^)f^ = —— - -— (r^ + r?)E^ 


(a^ — A)(aE — L)^ k+i 

+ ^- --2aE{L-aE)-A— =0. 

(r^ _l_ 

Now, as in the null geodesics, let us consider the special case L = aE and the general case separately. 

(a) Special case, L = aE-. 

Eor this case, the above equation becomes 

while Eqs. (A5a) and (A5b) become the same as for the null geodesics (A8a) and (A8b). 

(b) The general case [L — aE = x): 

Erom Eq. (AlO) with employing a reciprocal radius function (r = i) and substituting L — aE = x we hnd 

^‘^[2a;-l-fc(2a;+3)](j y2^2^2 ^-4^2 ^ ^ ^2^2^2 _ 2aExU^{l + U^U^) 

-M^(l + + (a^ - A)u^x^ = 0, 

where A = — '2^ + a?- + — ri^ for KNTN solution. The derivative of the above equation becomes 

271^4(2^2(1 + n^v^) - 2au^xE{l + 2n^u^) - + 2n^u^) + 2(a^ - A)u^x^, 

-(mu-1) ^ ^ ^ (1 + n 27 ( 2 )A *^^~^)/2 _|_ ^ 2^2 ^ 


Equations (A17) and (A18) can be combined to give 

(l+n2M2)2^2^y4^2(^2_^)_(^y_^)(^^^2y2) a('=+1)/2 (1 + 7727 ( 2 ) _ ^2^2 ^ 


and consequently 

2aEu^x(l + 7727 ( 2 ) = 2(a2 - A) 77 V - 7 ( 2(1 + 772772 )A(i +'=)/2 

- (77777 - 1)(1 + V ? U ^) —^(1 + 7727(2)A(''“^^/2 _)_ y2^2 

By eliminating E between these equations, we obtain the following quadratic equation for u^x^: 

A u'^x'^ + B 772^2 + C = 0, 

where we have defined 


+ (mu — 1)(1 + 772772 )]^ — 4 Aa 277 "^|, 

(A22) 

+ 772772 )! [ 2 A 772 + (77777 — 1)(1 + 772772 )]^ — 4 Aa 277 ^ 


- 1)(1 + 772772 ) [(77777 — 1)(1 + 772772 ) + 2 A 772 ] 1, 

(A23) 

2)2 [ 2 A 772 + (k + l)(mu — 1)(1 + 772772 )] ^ . 

(A24) 
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The discriminant V = — AAC of this equation is 

V = 64:A'^+^u^a^{l + + k‘^{mu - 1)2(1 + 712^2)2 

— [2Au2 + {mu — 1)(1 + v?u^)Y' I- 
The solution of the Eq. (A21) is 


(A25) 


2 2 

X U = — 


^(k 1)/2^2^ _|_ j^2y2^ 


^ k'{Z+ + Z-) ^ 2^/Aau^^Jk'‘^ - Z+Z_ 


2Z+Z_ 




(A26) 


where 

k' = k{mu — 1)(1 + n^u^), 

Z+Z_ = (2Au2 + {mu — 1)(1 + 'n?u’^)Y — AAa^u'^, 
Z± = 2Au^ + {mu — 1)(1 + v?u^) ± 2\/Aau‘^, 
and from the above equation we conclude that 


a; = ± < — 


/\(fc ^)/2(l -I-77,2^2^ 


^ k' {Z+ + Z-) ^ 2^/A au^yjk''^ - Z+Z_ 


2Z+Z_ 


Z+Z_ 


1/2 


(A27) 


When we put the values of x into Eq. (A20), we obtain the energy of circular orbit 

/ 

1 


E = 


— A 2 (fc + 1 )(to'u — 1)(1 + n2M2) 


^j2{l + rBu‘^) \ 

{mu — 1)(1 + ri^u^) — u^{a? — A) 




- Z+Z_ + k'{Z+ + Z_)/2 ± 2au^^^A^Jk'^ - Z+Z_ 


1/2 


(A28) 


In order to determine the ISCO of this spacetime, we need the explicit form of the effective potential. From Eq. 
(A15) we obtain 


1 -2 

-r = 


- 1 


+ ^e//> 


where 


^e// = 


1 


2' 2A-^P“-i-'=(2c^+3)] 

{a"^ - A){aE - Lf 2aA(L - aA) + A('=+i)/2 


2A^[2(.-;-l-fc(2i.c;+3)] [ 


1 + 


{A + rBy 


{A + V?) 


(A29) 


(A30) 


When the radial derivative of the effective potential vanishes, the particle stays in the circular orbit. For the ISCO 
equation, the second radial derivative of Veff must also vanish, i.e., = 0. We consider the region outside the 

“horizon” where A ^ 0. Then we have 


—A^A{n‘^ _l_ ^ 2 )| _ 2A2(n2 — 3r2) + A{k + 1) [n'^ + (4m — 3r)r] {n^ + r^) + (fc2 — l)(m — r)2(n2 + r2)2|. 

TdaEa; A3/2(n2 + r2)(n2 - 3r2) - 2x'^A^/'^ [2a^{n'^ - 5r2) - 2A(n2 - 5r2) + {n"^ + r‘^){n‘^ + 8mr - 7r^)] = 0.(A31) 
When we set Ex and x"^ using the expression above, we obtain the ISCO equation 


{k + l)(n2 + r2)2 [—Am — r{k — l)(m — r)2] + A[ — 4r^(a2 — A) — 4r^ — 4n2r^ + m{n^ + 5r‘^){n‘^ + r"^)] 
x| [2Ar + (m — r){n'^ + — AAr'^a^ + k{m — r){n'^ + r^) [2Ar + (m — r)(n2 + r2)] 


±2ar\fA 


Aa^r'^A — [2Ar + (m — r){r'^ + n^) ^ + k^{m — r)‘^{n^ + r^Y 


-11/2 


= 0 . 


(A32) 
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Having found the general expression, now let us consider a few special cases. First, we consider the radius of ISCO 
for the “extremal BD-KNTN spacetime,” where the parameters are related by 

(A33) 

in which the ’’horizons” of BD-KNTN spacetime 

r± = m± \Jrn? +v? — — a? (A34) 


coincides with the extremal Kerr radius r = m. Defining the reduced quantities 


1 

m 


n q 

— 5 Qm — ■> 

m m 


= — = \/l + n‘^-q'^, A = TO^(1 - r^)^, 
m 


we can express the extremal ISCO equation for BD-KNTN spacetime as follows: 

-{k + l)(r^ - 1)2 [l + {k- l)r^] {nl + rlf 

- I)"* [n^ -b 2n^r^(3 - dr^) -b r^(4g^ - 3r„)] 

X -b l)n))j -b 2n^rm(2 + k- dr^) -b Vm {rm [Mm + ?'m(-4 -b 2fc -b - fer^)] } 


(A35) 


±2rm\/l + nl, - ql. 


{k'^ - l)n^ -b 2n^rm [-2 + r^{k‘^ + 3)] -b {-4g^ -b r™ [4 -b rm(fc2 - 1)] } 


1/2 


= 0 . 


(A36) 


The solution of the above equation is Vm = 1 or r = m which coincides with the ISCO of extremal Kerr spacetime. 
For this special case, the ISCO radius of the BD-KNTN solution does not depend on the electrical and NUT charges 
or the Brans-Dicke scalar field. 

For another special case, we can calculate the ISCO equation for the BD Kerr solution by setting n = 0 and q = 0. 
Then we have 


{k + I) [—Am — r{k — l)(m — r)^] -b Ar( — da^ -b 4A — dr^ -b 5mr) 

X [(2A + r( m — r))2 — da^A -b k r{m — r) [2A — r(m — r)] ± 2aVA 

X =0, (A37) 

where here A = — 2mr -b . Our results are in accordance with corresponding solutions in GR [58-60]. 


[1] C. Brans and R. Dicke, Phys. Rev. 124, 925 (1961). 

[2] C. H. Brans, Phys. Rev. 125, 2194 (1962). 

[3] S. W. Hawking, Commun. Math. Phys. 25, 167 (1972). 

[4] M. Campanelli and C. O. Lousto, Int. J. Mod. Phys. D 2, 451 (1993). 

[5] K. A. Bronnikov, C. P. Constantinidis, R. L. Evangelista, and J. C. Fabris, Int. J. Mod. Phys. D 8, 481 (1999). 

[6] A. Bhadra and K. Sarkar, Gen. Relativ. Gravit. 37, 2189 (2005). 

[7] H. Kim, Phys. Rev. D 60, 024001 (1999). 

[8] D. H. Park, J. Korean Phys. Soc. 51, 258 (2007). 

[9] E. Barausse and T. P. Sotiriou, Phys. Rev. Lett. 101, 099001 (2008). 

[10] V. Cardoso, I. P. Carucci, P. Pani, and T. P. Sotiriou, Phys. Rev. D 88, 044056 (2013). 

[11] F. J. Ernst, Phys. Rev. 167, 1175 (1968). 

[12] F. J. Ernst, Phys. Rev. 168, 1415 (1968). 

[13] S. Chatterjee and S. Banerji, Acta Phys. Pol. B 11, 493 (1980). 

[14] T. Tsuchida and K. Watanabe, Prog. Theor. Phys. 101, 73 (1999). 

[15] M. Astorino, Phys. Rev. D 87, 084029 (2013). 

[16] H. A. Buchdahl, Phys. Rev. 115, 1325 (1959). 

[17] A. I. Janis, D. C. Robinson, and J. Winicour, Phys. Rev. 186, 1729 (1969). 

[18] H. A. Buchdahl, Int. J. Theor. Phys. 6, 407 (1972). 

[19] C. B. G. McIntosh, Commun. Math. Phys. 37, 335 (1974). 

[20] B. O. J. Tupper, Nuovo Cimento Soc. Ital. Fis. 19B, 135 (1974). 

[21] R. N. Tiwari and B. K. Nayak, Phys. Rev. D 14, 2502 (1976). 

[22] A. Eris and M. Gurses, J. Math. Phys. (N.Y.) 18, 1303, (1977). 



18 


[23] B. K. Nayak and R. N. Tiwari, J. Math. Phys. (N.Y.) 18, 289 (1977). 

[24] T. Singh and L. N. Rai, Gen. Relativ. Gravit. 11, 37 (1979). 

[25] H. Kim and H. M. Lee, Int. J. Mod. Phys. A 20, 6461 (2005). 

[26] A. Baykal and O. Delice. Gen Relativ. Gravit. 41, 267 (2009). 

[27] L. Akyar and O. Delice, Eur. Phys. J. Plus, 129, 226 (2014). 

[28] J. F. Plebanksi and M. Demianski, Ann. Phys. (N. Y.) 98, 98 (1976). 

[29] R. Debever, Bull. Soc. Math. Belg. 23, 360 (1971). 

[30] W. Kinnersley, J. Math. Phys. (N. Y.) 10, 1195 (1969). 

[31] C. Charmousis, D. Langlois, D. A. Steer, and R. Zegers, J High Energy Phys. (02) (2007) 064. 

[32] M. Astorino, J. High Energy Phys. 06 (2012) 086. 

[33] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, and E. Hertl, Exact Solutions of Einstein’s Field Equations, 2nd 
ed. Cambridge Monographs on Mathematical Physics (Cambridge University Press, Cambridge, England, 2003). 

[34] W. B. Bonnor, Z. Phys. 190, 444 (1966). 

[35] R. M. Misra and D. B. Pandey, J. Math. Phys. (N. Y.) 13, 1538 (1972). 

[36] H. Cebeci and T. Dereli, Int. J. Mod. Phys. D 13, 137 (2004). 

[37] F. J. Ernst, J. Math. Phys. (N. Y.) 17, 54 (1976). 

[38] F. Dowker, J. P. Gauntlet!, D. A. Kastor and J. Traschen, Phys. Rev. D 49, 2909 (1994). 

[39] F. J. Ernst and W. J. Wild, J. Math. Phys. (N. Y.) 17, 182 (1976). 

[40] G. W. Gibbons, A. H. Mujbata, and C. N. Pope, Classical Quantum Gravity 30, 125008 (2013). 

[41] R. Emparan, Phys. Rev. D 61, 104009 (2000). 

[42] C. Romero and A. Barros, Phys. Lett. A 173, 243 (1993). 

[43] N. Banerjee and S. Sen, Phys. Rev. D 56, 1334 (1997). 

[44] V. Faraoni, Phys. Rev. D 59, 084021 (1999). 

[45] A. Bhadra and K. K. Nandi, Phys. Rev. D 64, 087501 (2001). 

[46] D. K. Ciftci and O. Delice, J. Math. Phys. (N. Y.) 56, 072502 (2015). 

[47] P. Kirezli, D. K. Ciftci, and O. Delice, Gen. Relativ. Gravit. 45, 2251 (2013). 

[48] B. Chauvineau, Classical Quantum Grav. 20, 2617 (2003). 

[49] B. Chauvineau, Gen. Relativ. Gravit. 39, 297 (2007). 

[50] N. N. Bocharova, K. A. Boronnikov, and V. N. Melnikov, Vestn. Mosk. Univ. Fiz. Astron. 6, 706 (1970). 

[51] J. D. Bekenstein, Ann. Phys. (N. Y.) 82, 535 (1974). 

[52] J. D. Bekenstein, Ann. Phys. (N. Y.) 91, 75 (1975). 

[53] A. Anabalon and H. Maeda, Phys. Rev. D 81, 041501 (2010). 

[54] A. Anabalon, J. High Energy Phys. 06 (2012) 127. 

[55] S. Bhattacharya and H. Maeda, Phys. Rev. D 89, 087501 (2014). 

[56] Y. Bardoux, M. M. Caldarelli, and C. Charmousis, J. High Energy Phyds 05 (2014) 039. 

[57] M. Astorino, Phys. Rev. D 91, 064066 (2015). 

[58] C. Chakraborty, Eur. Phys. J. C 74, 2759 (2014). 

[59] P. Pradhan, Classical Quantum Gravity 32, 165001 (2015). 

[60] S. Chandrasekhar, The Mathematical Theory of Black Holes, International Series of Monographs on Physics Vol. 69 (Oxford 
University, New York, 1992). 

[61] C.-M. Claudel, K. S. Virbhadra, and G. F. R. Ellis, J. Math. Phys. (N. Y.) 42, 818 (2001). 



